Eigenstates of Operating Quantum Computer: 
Hypersensitivity to Static Imperfections 
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We study the properties of eigenstates of an operating quantum computer which simulates the 
dynamical evolution in the regime of quantum chaos. Even if the quantum algorithm is polynomial 
in number of qubits n q , it is shown that the ideal eigenstates become mixed and strongly modified 
by static imperfections above a certain threshold which drops exponentially with n q . Above this 
threshold the quantum eigenstate entropy grows linearly with n q but the computation remains 
reliable during a time scale which is polynomial in the imperfection strength and in n q . 
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Feynman suggested that a quantum computer could 
simulate quantum mechanical systems exponentially 
faster than a classical computer |ÏJ while Shor signif- 
icantly extended this class by his ground-breaking al- 
gorithm for integer factorization More recently, a 
few quantum algorithms which achieve the exponential 
speedup have been developed for various quantum and 
classical physical systems, ranging from some many-body 
problems || to spin lattices 0], and models of quan- 
tum and classical chaos These algorithms have 
been constructed for idealized perfect quantum comput- 
ers. In reality, any laboratory implementation of a quan- 
tum computer would involve imperfections B. The first 
investigations have shown a certain stability of quantum 
evolution and algorithms with respect to decoherence ef- 
fects JqJ , noisy gates JÏ0|-^2|, and static imperfections 
|j],U ■ These studies have focused on the fidelity of quan- 
tum computation as a function of time during the real- 
ization of a given quantum algorithm. 

In this Letter, we study the properties of the eigen- 
states of an operating quantum computer in the presence 
of static imperfections. The computer is simulating effi- 
ciently the time evolution of a dynamical quantum sys- 
tem described by the sawtooth map [Q . We focus on the 
regime of quantum ergodicity, in which eigenfunctions 
are given by a complex superposition of a large number 
of quantum register states. In this regime, the effect of 
a perturbation is enhanced by a factor which is exponen- 
tial in the number of qubits. This phenomenon has close 
links with the enormous enhancement of weak interac- 
tions in heavy nuclei In the following we illustrate 
this general effect for the case of static imperfections in 
a realistic model of quantum computer hardware. 

The classical sawtooth map is given by 

n = n + k{6 -n), 6 = e + Tn, (1) 
where (n,8) are conjugated action-angle variables (0 < 



< 27r), and the bars denote the variables after one 
map step. Introducing the rescaled momentum variable 
p = Tn, one can see that the classical dynamics depends 
only on the single parameter K = kT, so that the mo- 
tion is stable for —4 < K < and completely chaotic for 
K < —4 and K > 0. The quantum evolution for one map 
iteration is described by a unitary operator U acting on 
the wave function tp: 



-iTn^/2 ik(0—Ky/2 



(2) 



where n — —íd/d0 (we set fi = 1). The classical limit cor- 
responds to k — > oo, T — > 0, and K — kT — const. In this 
Letter, we study the quantum sawtooth map (||) in the 
regime of quantum ergodicity, with K = —n<p<TT 
(torus geometry). The classical limit is obtained by in- 
creasing the number of qubits n q — log 2 N (N number of 
levels), with T = 2n/N (k = K/T, -N/2 <n< N/2). 
The quantum algorithm simulates with exponential 
efficiency the quantum dynamics (|^) using a register of 
n q qubits. It is based on the forward/backward quantum 
Fourier transform |Ï4| ] between the 6 and n representa- 
tions and requires 2n q Hadamard gates and 3n q — n q 
controllcd-phase-shift gates per map iteration |?J. 

Following we model the quantum computer hard- 
ware as an one-dimensional array of qubits (spin halves) 
with static imperfections, i.e. fiuctuations in the individ- 
ual qubit energies and residual short-range inter-qubit 
couplings . The model is described by the many-body 
Hamiltonian 



H 8 = J2(&o + 5 i )àl 



X! ■'■j (j · (j j- 

i<j 



(3) 



where the &i are the Pauli matrices for the qubit i, and 
A is the average level spacing for one qubit. The second 
sum in (j^) runs over nearest-neighbor qubit pairs, and 
8i , Jij are randomly and uniformly distributed in the in- 
tervals [—5/2,6/2] and [— J, J], respectively. We study 
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numerically the many-qubit eigenstates of the quantum 
computer (JsJ) running the quantum algorithm described 
above. The algorithm is realized by a sequence of instan- 
taneous and perfect one- and two-qubit gates, separated 
by a time interval r g , during which the Hamiltonian (j^) 
gives unwanted phase rotations and qubit couplings. We 
assume that the average phase accumulation given by Ao 
can be eliminated, e.g. by means of refocusing techniqucs 

0; 

Since the evolution operator (H) remains periòdic in the 
presence of static imperfections, Ü^(t + T) = U^(t) 
(e = 6r g ), all the informations about the system dynam- 
ics are included in the quasienergy eigenvalues Aq ' and 
eigenstates <pa of the Floquet operator: 

Í>«(T)^)=exp(^)^. (4) 

In Fig| (top left) we show the parametric depcndcncc 
of the quasienergy eigenvalues on the dimensionless im- 
perfection strength e, for a given realization of at 
J = and for n q — 9 qubits. One can clearly see the pres- 
ence of avoided crossings, a typical signature of ergodic 
dynamics. The variation of a given quasienergy eigen- 
state with e is illustrated by the Husimi functions |Ï(J of 
Fig.jï]. At e = the eigenfunctions display a complex pat- 
tern delocalized in the phase space (see Fig.^top right). 
The symmetries of the Husimi functions (8 — ► 2tt — 9, 
p — ► — p) are destroyed when e ^ |Ï7| . The eigenfunc- 
tions in the presence of imperfections give a good repre- 
sentation of the unperturbed (e = 0) eigenfunctions of the 
quantum sawtooth map (g) at most until the first avoided 
crossing. After that one cannot make a one to one corre- 
spondence with the unperturbed case. This is confirmed 
by the last two Husimi functions of Fig.|], taken for the 
chosen level in the vicinity of the first avoided crossing 
(bottom left, e = 4 x 10~ 4 ) and for a stronger imper- 
fection strength (bottom right, e = 10~ 3 ). In the first 
case there is still some similarity with the corresponding 
Husimi function at e — 0, whilc in the latter case any 
resemblance has been effaced. 

A more quantitative indication of the similarity be- 
tween exact and perturbed eigenstates is provided by the 
quantum eigenstate entropy, 

JV 

S a = - ^2 Vad log 2 Pa/3, (5) 
(3=1 

where p a p = |(</^ 0) |</>i £) )| 2 . In this way S a = if <fi^ 

coincides with one eigenstate at e = 0, S a = 1 if 4$ 
is equally composed of two ideal (e = 0) eigenstates, 
and the maximal value S a — n q is obtained if all </>L 

((3 = 1,...,N = 2 n ") contribute equally to (f>$ . In or- 
der to reduce statistical fluctuations, we average S a over 
a = 1, ...,N and over 3 < Np < 10 3 random realizations 
of Si , Jij . In this way the total number of eigenstates is 



NpN rí 10 4 . The variation of the average quantum en- 
tropy S with e is shown in Fig.S. It demonstrates that S 
grows from S' = 0ate = 0toa saturation value S ~ n q 
corresponding to maximal mixing of unperturbed eigen- 
states. We study this crossover for 4 < n q < 12 at J = 
and find that the mixing takes place at smaller vàlues 
when n q increases. In Figj^ we show that the critical 
imperfection strength e x at which S = 1 drops expo- 
nentially with the number of qubits. This exponential 
dependence holds also in a simple toy model with a sin- 
gle impurity, in which an energy fluctuation S is switched 
on for a single qubit and only for one time interval r g be- 
tween two elementary gates (e.g., after the first quantum 
Fourier transform). 




FIG. 1. Parametric dependence of the quasienergy eigen- 
values on the imperfection strength e for a given random re- 
alization of Si, at J = 0, n q — 9 (top left); for the red-colored 
level the corresponding Husimi functions in action-angle vari- 
ables (p, 9) (— 7r < p < n -vertical axis- and < 9 < 2tt -hori- 
zontal axis-) are given at e = (top right), e = 4 x 10 -4 (bot- 
tom left), and e = 10~ 3 (bottom right). We choose the ratio 
of the action-angle uncertainties s — Ap/A9 — TAn/ A9 — 1 
(ApA9 = T/2). The color is proportional to the density: blue 
for zero and red for maximal density. 

The exponential drop of the threshold e x can be un- 
dcrstood following a theory originally developed for the 
parity breaking induced by weak interaction in the scat- 
tering of polarized neutrons on complex nuclei [Ï3]| . In- 
deed, due to quantum chaos in (||), the Floquet problem 
(|]) has ergodic eigenstates, (j)a ^ = ELi^" 1 '' 1 "" wrtn 
u m being quantum register states and Ca randomly fluc- 
tuating components, with |ci" l ' l | ~ l/\/N. For the model 
with a single imperfection Sàf , acting on a time interval 
T g , the transition matrix elements between unperturbed 
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eigenstates have a typical value: 

N 
m— 1 

(6) 

The last estimate for V4 yp results from the sum of N un- 
correlated terms. Since the spacing between quasienergy 
eigenstates is Ai? ~ 1/iV, the threshold for the breaking 
of perturbation theory can be estimated as 

V typ /AE ~ e x VN ~ 1. (7) 

The analytical result ~ 1/ \HV is confirmed by the 
numerical data in Fig.g. The same theoretical argument 
gives an exponential drop of e x for the static imperfec- 
tion model (||). In this case, the estimate can be ob- 
tained with S — > (sum of n q random detunings Si) 

and T g — ► T g n g ~ Tg-ri^. This gives e x ~ N~ l / 2 n q J ^ 2 , 
again in good agreement with the data of Fig]|. For the 
case J = <5, the threshold e x decreases by a factor sa 1.5 
at n q — 9 with respect to the J = case (see Fig.|]), 
since additional qubit couplings are introduced. We note 
that the hypersensitivity to perturbations has been pro- 
posed as a distinctive feature of chaotic dynamics [fill. 
However, the authors of Ref. [Q considered the effect 
of a stochastic environment, while we consider a closed 
Hamiltonian system. 
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FIG. 2. Quantum eigenstate entropy 5* as a function of the 
scaled imperfection strength e for J = 0, n q — 5 (triangles), 
7 (diamonds), 9 (empty squares), and 11 (circles), for J — 5, 
n q — 9 (filled squares), and for the single imperfection model 
at n q = 11 (stars). The straight lines give the theoretical 
estimates 2 S = Ae 2 N (sòlid line) and 2 S = Be 2 n 5 q N (taken 
at n q = 11, dashed line), with the numerically determined 
constants A = 0.37 and B — 0.25. 

In the mixing regime (e > e x ) the number M of unper- 

turbed eigenfunctions 4>^ , which have a significant pro- 

jection over a given <pa\ is exponentially large. For the 



single imperfection model one has M ~ 2 S ~ T/AE ~ 
e 2 N, since the mixing takes place inside a Breit-Wigner 
width given by the Fermi golden rule: T ~ V^ p /AE ~ e 2 . 
The above estimate for M is in agreement with the nu- 
merical data of Fig.|| and is similar to that one used in 
HPl for onset of quantum chaos in the static model (||). 
We emphasize that this estimate implies that the quan- 
tum eigenstate entropy grows linearly with the number 
of qubits n q . In the Fermi golden rule regime, the life- 
time (measured in number of kicks) of an unperturbed 
eigenfunction is given by r x ~ 1/r oc 1/e 2 p||l·1| . If the 
imperfections are described by the model ffi), one has 
t x r~j l/(e 2 n^). Therefore a reliable quantum computing 
of the dynamical evolution of the model (U) is possible 
up to a time scale which drops only algebraically with 
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FIG. 3. Dependence of the coupling e x at which S = 1 
on the number of qubits, for J — (circles) and for 
the single imperfection model (squares). The lines give 
the theoretical dependences e x = A -1 ^ 2 N -1 ^ 2 (above) and 
e x — B^ 1 / 2 iV -1 / 2 ^ 5 / 2 (below), with the constants A and B 
obtained from the data of Fig.0. 

the number of qubits, in agreement with the findings of 
Ref. JtJ. In FigJJ we show the fidelity of quantum evo- 
lution, f(t) = |(^(°)(í)|^(í))| 2 . In the top figures, the 
initial state is an unperturbed eigenstate. For e < e x 
(S < 1), the fidelity is very close to 1 at all times, 
since the eigenstates are not mixed by the imperfections 
(see Fig.[|(a)). On the contrary, in the Fermi golden 
rule regime e >> e x a perturbed eigenstate, when de- 
composed into the unperturbed eigenstates, contains a 
large number of components (S 3> 1). The distribution 
of these components over energy, called local density of 
states, has a typical Breit-Wigner shape of width T. Since 
its Fourier transform drives the fidelity decay p||Ï9|] , one 
obtains f(t) w exp(— Tt), in agreement with the data of 
Fig.||(b) . The exponential decay continues up to a value 
/ w Í/2 S given by the inverse of the number of levels 
mixed inside the Breit-Wigner width. The case in which 
the initial wave function ip(0) is a momentum eigenstate 
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is considered in Fig.[|(c,d). In this case ip(0) projects 
significantly over order N unperturbed eigenfunctions. 
Therefore, for e < e x , f(t) displays a Gaussian decay (see 
Fig.^(c), and also Ref. 0): in this regime the imperfec- 
tions do not change significantly the eigenstates (S < 1), 
but the initial state is composed of many eigenstates and 
a Gaussian decay of f(t) is espected from perturbation 
theory 20|2ÏJ]. FigJ^(d) shows that in the Fermi golden 
rule regime (e > e x ) the fidelity decays exponcntially, 
with rate T given by the Breit-Wigner width ]2Ï] ]. The 
decays stops when / ss Í/N, namely when / approaches 
the inverse of the dimension of the Hilbert space. 
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FIG. 4. Fidelity as a function of time, for n q — 9 qubits, 
J = 0, e = 10~ 4 (left) and e = 3 x 10~ 3 (right), with ini- 
tial wave function a Floquet eigenstate at e = (top) or a 
momentum eigenstate (bottom). The dashed lines show the 
exponential decay f(t) — exp(— t/tf), with tf ~ 7. 

In summary, we have shown that the eigenstates of a 
quantum computer simulating a system with quantum 
chaos are hypersensitive to static imperfections: they 
are significantly different from the exact eigenfunctions 
above an imperfection strength threshold which drops 
exponentially with the number of qubits. This border is 
directly relevant for quantum algorithms which aim at 
computing static properties of physical systems, for ex- 

We also 
assume the 



ample energy eigenvalues and eigenvectors 2 
note that quantum adiabatic algorithms |23 
ability to drive the evolution of a quantum state towards 
the ground state of a specified Hamiltonian. Moreover, 
a few relevant quantum algorithms can be formulated 
in terms of determining the eigenvalues corresponding 
to eigenstates of a given unitary operator [^.^5). Fur- 
ther investigations are required to establish under which 
conditions the observed hypersensitivity of eigenstates to 
perturbations could be relevant for the stability of these 
algorithms. 
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